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Abstract 

Dynamics near the grazing manifold and basins of attraction for a 
motion of a material point in a gravitational field, colliding with a mov- 
ing motion-limiting stop, are investigated. The Poincare map, describing 
evolution from an impact to the next impact, is derived. Periodic points 
are found and their stability is determined. The grazing manifold is com- 
puted and dynamics is approximated in its vicinity. It is shown that on the 
grazing manifold there are trapping as well as forbidden regions. Finally, 
basins of attraction are studied. 
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1 Introduction 

Vibro-impacting systems are interesting examples of non-linear dynamical sys- 
tems, exhibiting new levels of complicated dynamics due to their non-smoothness 
and have important industrial applications, see [1] and references therein. Re- 
search in this field was initiated by studies of impacting linear oscillators [5] - 
[20] as well as closely related accelerator models of particle physics originating 
from the Fermi model [5T] - [35] . A very characteristic feature of such systems is 
the presence of non-standard bifurcations such as border-collisions [2] - [S] and 
grazing impacts [8], [14] - [20] . Such events often lead to complex chaotic motion, 
see [T] for a review of border-collision and grazing bifurcation phenomena. Little 
work was done, however, on basins of attraction of different attractors typically 
coexisting in impacting systems. Early results were obtained in [141 1151 [T7] . It 
seems however, that the structure of basins of attraction, especially in the case 
of coexistence of grazing dynamics and periodic orbits, deserves a further study. 

In the present paper we investigate dynamics near the grazing manifold 
and basins of attraction for a motion of a material point in a gravitational 
field colliding with a moving motion-limiting stop. Typical example of such 



dynamical system, related to the Fermi model, is a small ball bouncing vertically 
on a vibrating table [29] - [32]. Since evolution between impacts is expressed 
by a very simple formula the motion in this system is easier to analyze than 
dynamics of impact oscillators. 

The paper is organized as follows. In the second Section of this paper a one- 
dimensional motion of a material point in a gravitational field, colliding with a 
limiter, representing unilateral constraints, is considered. It is assumed that, at 
impact, there is an inelastic rebound with coefficient of restitution < R < 1 
[33] [33]. Next, the Poincare section is defined and the return map, describing 
evolution from an impact to the next impact, is derived generalizing and extend- 
ing our previous results [35] . In Section 3 we find periodic points and determine 
their stability. In Section 4 the grazing manifold is computed. It turns out that 
the Poincare map is not well defined on the grazing manifold. However, we 
demonstrate that the map can be continued to the whole grazing manifold after 
appropriate regularization procedure, except from some critical points. Then 
the nonlinear equation for time of the next impact is approximately solved and 
dynamics is linearized in the vicinity of the grazing manifold. It is shown that 
on the grazing manifold there are trapping as well as forbidden regions which 
are determined explicitly. Finally, basins of attraction are studied in Section 5. 
More exactly, we investigate boundary between basins of grazing dynamics and 
periodic motion which, for large velocities, becomes very complicated. 

2 Motion with impacts: the Poincare section 
and return map 

Let us consider motion of a material point moving vertically in a gravitational 
field and colliding with a moving motion-limiting stop, representing unilateral 
constraints. Let motion of the point in time intervals between impacts be de- 
scribed by equation: 

mx = —mg, (1) 
where x = 4|, while motion of the limiter is expressed as: 

y = y(t), (2) 

with known function y(t). We assume that y, y, y are continuous functions of 
t for t G [0, oo). 

Equations of impact are of form: 

x{Ti)=y{Ti), (3a) 
i (rf) - y (r<) = -R [x (rr) - y ( Ti )] , (3b) 

where it is assumed that the limiter's mass is infinite. In ([3]) Ti denotes the in- 
stant of i-th impact, t~ (t^~) is the time just before (after) the i-th. impact and 
R is the coefficient of restitution, < R < I [351 134] . The coefficient of resti- 
tution is a measure of the elasticity of the collision between the material point 
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and the limiter, and is also a measure of energy dissipated during a collision. 
The coefficient of restitution can be determined experimentally or computed on 
the basis of a specific model of collisions. In the present work we shall assume 
that R is known and constant. 

It follows from the dynamics that x(t) is a continuous function of time while 
x (t), continuous within all time intervals (t^Lj, t~) , is subject to discontinuous 
changes in time instants Tj described by the impact equation (|3b[) . By x (t^~), 
x (rf) we understand the left-sided and right-sided limits of x (t) for t — > Ti, 
respectively, for every i = 0, 1, . . . . In what follows we shall not assume a priori 
knowledge of Tf's. 

We shall now consider equation (JTJ for t G { T t ■> rjlj). General solution 
reads: 

x (t) = -l 9 t 2 + 4h + 4\ (4) 

x(t) = -gt + cf, (5) 

where c{ , are arbitrary constants. 

The constants cj^, in (gj), © can be determined from the initial condi- 
tions x(t+), x(t+): 

cf = x{ri)+gT h c« = y^-c^r l +x(r l ). (6) 
The solution (|U), ([5|) can be now rewritten in form: 

x(t) = ~\gt 2 + [x (r+) + gn] t - \grj - x (r+) r, + x(r,), (7a) 
x(t) = -gt + x(r+) +gu, (7b) 

for each i = 0, 1, ...and £ 6 (t+, t~ | _ 1 ). Substituting £ = r i+1 in (I7a|) we get the 
formula for the position of the material point after the i + 1-th impact: 

x (t 1+ i) = x (Ti) - \g (r i+1 - t^ 2 + x (r+) (r i+1 - Tj) . (8) 

It is now possible, applying to §5§ the impact condition (|3ap. to obtain a 
nonlinear equation from which the time of the next impact can be computed 
from the given impact time and the corresponding velocity: 

V (Ti+i) = V (Ti) - \g (r i+ i - r,) 2 + x (r+) (r i+ i - r») . (9) 

On the other hand, substituting t = t^ +1 in (|7bp we obtain the formula for 
the velocity of the material point just before i + 1 -th impact: 

* ( T i+i) = * ( T 4 + ) - 9 {ri+i - Ti) . (10) 
To obtain the velocity just after the i + 1-th impact the equation (|3b[) is used: 

x (r+ = -itt (r H1 ) + (1 + R) y (r i+1 ) . (11) 
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Combining equations (fTTj) and (|10[) we obtain: 

x (r+ x ) - i? [-i (r+) + .g (r i+ i - r<)] + (1 + R) V (r i+ i) . (12) 

Finally, equations mapping the time of the i-th impact to the time of the 
i + 1-th impact as well as the velocity of the material point just after the i-th 
impact to the velocity just after the i + 1-th impact for arbitrary motion of the 
limiter consist of Eqs.([U), [fT2]): 

y (n+i) = y (n) - \g (n +1 - t 4 ) 2 + ± (t+) (t. 1+1 - n) , (13a) 

* fa+i) = ( T t) + R 9 - n) + (1 + R)y (r i+1 ) . (13b) 

With help of Eqs. (fl3| and Eq.® the solution can be continued to the 
interval (rf +v tJ +2 ).^ 

Introducing non-dimensional variables: 

Ti = un, Y (T() = y (n) /a, ^ = fx (r+) , A = (14) 

where a; and a determine time and length scales, we obtain from (|13|) : 

AF (T m ) = AF (^) - (T i+1 - T,) 2 + Vi (T t+1 - Ti) , (15a) 
V l+l = -RV- + 2R (T t+1 - + A (1 + R) Y (T i+1 ) , (15b) 

where Y = 4|? . It follows that the map (fl5|) depends on two control parameters 
only: R and A. 

Assuming finally position of the limiter in form y (t) = a sin (art + ip), that 
is y (T) = sin (T + ip), we get: 

A sin (T i+1 + ip) = A sin (T< + - (T i+1 — Ti) 2 + V t (T l+1 - T<) , (16a) 
y !+ i = -JZV5 + 2i? (T m - Ti) + A (1 + i?) cos (T m + 99) . (16b) 

The map (| is invariant under the translation Tj — * 7$ + 27r and thus the 
phase space (T, V) is topologically equivalent to the cylinder. Accordingly, Eqs. 

(|16p take a simpler form if a new variable Xi = Ti + ip (mod 2n) is introduced: 

A sin (AQ +1 ) = A sin (X t ) - (AQ +1 - X,) 2 + V t {X i+1 - X,) , (17a) 
V i+1 = -RVi + 2R (X i+1 - Xi) + A (1 + R) cos (X i+1 ) . (17b) 

The map (|17p must meet several physical conditions to correspond exactly 
to the original physical problem described by equations (fTJ) , @, ([3]). First of 
all, we have to choose the right solution out of infinitely many solutions of the 
nonlinear equation (|17a|) . We thus choose the solution X i+1 with the smallest 
nonnegative difference Xi + i — JQ. Moreover, since Vi is the velocity of the 
material point just after the impact, it must not be smaller than the velocity of 
the limiter. Since the quantity 

^-Acos(A- i ) = fe(i(r+)-y(r i )) > (18) 
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is the relative velocity of the material point in non-dimensional units we imple- 
ment equations (|17[) with the following physical conditions: 

6i = X l+1 - X, > , (19a) 
Wi = V i -Xcos(X i ) > , (19b) 

which define the physical phase space (X, V). 

Equations similar to (flU)) were derived in for position of the limiter 
in the special form y (r) = a sin (lot + ip) and correct stability conditions for 
fixed points of the map were determined. However, the authors did not use 
the non-dimensional variables as in (|15[) or (|17p and hence their map contained 
redundant control parameters. Moreover, they did not compute the map for 
the velocity just after the impact, x (t^~), for which the important inequality 
x (rf) > y (rj) or (|19b[) holds, but for the velocity just before the impact, 

Computations of impact times TVs from nonlinear equation (|16a[) are com- 
plicated. The problem simplifies greatly for motions for which the approximate 
equality x (Ti + i) = x(Ti) or, equivalently, y(T J+ i) = y (Ti) holds. The condi- 
tion y (Ti + i) — y (Ti) leads to A sin (Ti + i + ip) = A sin (Ti + ip) which applied to 
(|16ap yields Ti + i — Ti = Vi. Hence Eqs.fTTTJ]) reduce to the following approximate 
map: 

T t+1 =Ti + Vi, (20a) 
V-+i = RVi + A (1 + R) cos (Ti + Vi + <p), (20b) 

equivalent for ip = it to the map introduced by Holmes 29J (note that A (1 + R) 
is the control parameter 7 of Holmes) . The map (|2T)|) was studied theoretically 
[2"9"] while experimental results were reported in [3D1 [31] ■ 



3 Fixed points and their stability 

Since the map (fT?} is 2ir periodic in X the fixed points (1 - cycles) are found 
by substituting X^™\ — x\ n ^ = 2mr. Accordingly, we get: 

T4 (n) = 27m, (21a) 
cos(xi")) = ^# ) (21b) 

where LX* n \ V*"M is a fixed point on the cylinder S x K. In (|2T|l n is an 
arbitrary non-zero integer such that: 

-i^^irw 1 ^ 1 (" = ±1 ' ±2 ' ■■■)• ( 22 ) 

The periodic motion of the material point is thus possible when the condi- 
tions fl2T]), are fulfilled. 
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From (|21a[) and (|21b[) we have, for a given n, the pair of fixed points: 

Vi n) = 2™, X<$ = arccos (^f) , X^ = 2tt - X^ . (23) 



Let us notice here that the fixed points of the map (115|) are exactly the same 
as the fixed points of the approximate map (I20[) but, of course, their stabilities 
should differ. 

Stability of periodic solutions l|2"3"|) can be investigated by substituting in 
(fT7|) Xi = xi'p + di, Vi = V* + Uj with small deviations di, m from the fixed 
point. In all formulae below we have p = 1,2. After some rearrangements we 
get equations for the perturbations di, u,: 

A sin (d i+ i + Xl^) = A sin (di + X^ - (d l+1 - d t f - 4tto {d l+1 - d t ) + (24a) 

+ 2irn (di+i — eij) + 2irnUi, 
u i+1 = -Rut + 2R (d i+1 - di) + (24b) 

+ A (1 + R) [cos + - cos (x$A 



The system of equations (|24|) has the trivial solution d; = Ui = . Lineariz- 
ing Eqs. lfM)) around the trivial solution we obtain: 

Xd i+ i cos (xip^J — Xdi cos (xi^J — 2Tm(d i+ i — di) + 2-nnui, (25) 
u i+ i = -Rui + 2R{d l+1 - di) - A(l + R)d i+ i sin (xty) , (26) 



and substituting Xip from Eq. (|23[) we obtain the linearized map defined via 
the stability (Jacobi) matrix §: 

d i+1 \ =§ fdi\ 

Ui+1 J \ Ui J 

( 1 1±R \ 

§ = , , . ( 27b ) 



where /4™ 2 = iy'A^ (1 + R) - 4ir 2 n 2 (1 - R) . Let us note that due to d22j) 
(«) ' i 

/X^ 2 are rea l- 

The trivial solution di = Ui = is asymptotically stable if and only if the 
roots Ai, A2 of the characteristic equation of the stability matrix: 

f-A i±S \ 

det I , , / , = A 2 +a„A + (3 = 0, (28a) 

' -/!<»> (i?2_i±« M («))_ A j 

i±^4 n) -l-i? 2 , /3 = i? 2 , (28b) 
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fulfill inequalities |Ai i2 < 1. Using the Schur-Cohn criterion we get: 

[3<l, \a p \ </3+l. (29) 

Let us first consider the first fixed point , vj n ^ = Pi""*. For p = 1 

we obtain from Eqs.([2"8|). (|2"9")l conditions guaranteeing asymptotic stability of 
pin). 

R < 1 , (30a) 

A („) ^ 2jn£-R) < A < A („) ^ I _2j JI ^2 (1 _ fl 2 ) 2 + 4(1+iJ 2 ) 2 , (30b ) 

assuming that the inequality (f22|) is fulfilled. It turns out that stability of the 
fixed points of the exact map (|15[) is indeed different than the stability of these 
points for the approximate map (|20|) , see [29] . 

On the other hand, the second fixed point (X%>, K (n) ) =Pi n) is unstable. 

Indeed, for A > x£) and < R < 1 the discriminant of the characteristic 
equation (|2"5|) is positive so that Ai, A 2 are real. Furthermore, the inequalities 
< AxA 2 = R 2 , Ax + A 2 > 1 + R 2 result from (|2SJ) for p = 2 and A > A^. It 

follows finally, that < Ai < R 2 < 1 and 1 < A 2 . The fixed point is thus 
a saddle. This result and the inequality (|30b[) were also derived in [3"2"] . 

4 Grazing dynamics 
4.1 Bouncing ball 

The map (fTTf has also the manifold of fixed points, defined by the condition: 

X l+1 =Xi=X„ V l+1 =Vi=V* = XcosX*, (31) 

what is verified by the direct substitution into Eq. lfTT]) . We shall refer to the 
manifold (|3ip as the grazing manifold. In this Section we shall analyze dynamics 
near the grazing manifold. 

Let us first notice that a function Xi + i = f (Xj, Vj) is defined implicitly by 
Eq.|l7a|: 

F{X i+1 ,X h Vi) i AsinX l+1 -\smXi + (X i+1 - X t ) 2 + 
-(X i+1 -Xi)Vi = Q. 

Therefore since 

dF(X i+1 , X i} V t ) 



i+1 



A cos Xi+i + 2 (X i+1 - Xi) - Vi, (33) 



then on the grazing manifold dF (X i+ i, Xi, Vi) /dX i+ i = 0, see Eg. ([31) . and 
thus Xi + i — f (Xi, Vi) cannot be defined therein as follows from the implicit 
function theorem 1371. 
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The map can be however continued on the whole grazing manifold except 
from some critical points. Let us first notice that since there was impact at time 
Xi, it follows that X i+1 = Xi is a trivial solution of equation (j3"2")) , corresponding 
to the imposed initial condition. This trivial solution is the source of singular 
behaviour of the map on the grazing manifold. We can get rid of the trivial 
solution dividing the function F {X i+ i, Xi, Vi) by X i+i — Xi and expressing as 
a function of Vi, Xi and Si = — Xf. 



G(S U X, V^g (34a) 
G (5 U X h Vi) = j- (sin (5, cosX, + (cosS t - 1) sinX;) + S l -V t = 0. (34b) 

Let us notice that the function G (Si, Xi, Vi) is well defined, the limit- 
ing case Si — > including. Equation (|34[) defines implicitly a function Si — 
g(Xi, Vi). Since 

dG{5l d ^ Vl) = HAsinX, + 1) + O (Si) , (35) 

it follows that Si = g (Xi, Vi) can be defined everywhere on the grazing mani- 
fold, i.e. for Si — * 0, except from critical points: 

1 - ±\smX cr = 0, (36) 

which appear for A > 2, where the assumption of the implicit function theorem, 
dG(Si, X t , Vi) /dSi ^ 0, is not fulfilled. It follows easily from Eqs.([3l|), l[36|) 



that in critical points V cr (A) = ±v A 2 — 4. 

It is important that while the quantity Si is non- negative, see (|19a[) . it is 
also small near the grazing manifold: 

0<5i = X i+1 -Xi<l. (37) 

It is thus natural to expand G (Si, Xi, V.) in a power series of Sf. 

G(S h Xi, V) = (XcosX, - Vi) + (1 - iAsinJSQ) S % - ±AcosJQ5 2 -' 
+^AsinJQ5f + ... =0, 



(38) 



which converges under assumption (|37p . Very close to the grazing manifold 
Si <C 1 and we can neglect higher-order terms in l|38j) to get approximate formula 
for the function g (Xi, Vi): 

Si = g (X t , Vi) = ^\ cosX - + O (S*) . (39) 



Eq. (|3"9"|) in the initial variables (TH)) has the following form 

-2 



Ti+i - Ti = Vi=x (r+) - y (n) , Ai=x (r+) - y (n) , (40) 
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where Vi and At are relative post-impact velocity and acceleration, respectively. 
We shall demonstrate in the next Subsection that the formula (l40|) is more 
general. 

The approximation (|39p breaks down in critical points on the grazing mani- 
fold defined by (p?l)j) . that is v; = 0, Ai = 0. On the other hand, we can use (|3T))) 
near the grazing manifold wherever the consistency conditions 

< Si = ' | A "'- V < 1, (41) 

are fulfilled, see (|37|) . It follows from inequality (|19b| that the condition (|4lT) is 
violated for 

l-±AsinXi<0. (42) 

Therefore the region (|42p of the grazing manifold is non-physical and thus for- 
bidden, see Fig. 1. On the other hand, Eq. (|40]) shows that in the forbidden 
region Vi — and Ai > 0, c.f. the definitions (|40|) . and thus this region is 
repelling. The acceleration function was also used to analyze motion of impact 
oscillator [TH [HI US] . 

The map (|17[) can be simplified now in the neighbourhood of the grazing 
manifold. Equation (|17a[) can be approximated by ([39]) and we rewrite (fT7|) in 
the following form: 

Wi = \1 — ^AsmXiJ Si , (43a) 

w i+ i = —Rwi + 2R5i + XR (cos X i+1 — cos Xj) , (43b) 

where Wi = Vi — AcosX;, Si = X J+ i — Xi. Expanding now cos(X + i) = 
cos {Xi + Si) and neglecting higher-order terms in Si we obtain the map: 

Wi = (l - ^AsmX^ S l , (44a) 

w i+ i = -Rwi + 2R(l - |AsinXi) 5 U (44b) 

which, finally, takes a very simple form: 

w l+1 = Rw t . (45) 

Hence for R < 1 this part of the grazing manifold where the condition (|4ip 
is fulfilled is attracting. Therefore Wi — > and it follows from (|44ap that also 
Si -» for 1 - ±A sinXi — ► 1 — ^AsinX* > 0, i.e outside critical points. In the 
final stage of grazing collisions are thus so rapid that the motion of the limitcr 
can be neglected. 

We shall refer to the region on the grazing manifold 1 — ^AsinX < as the 
forbidden or repelling region while the attracting part of the grazing manifold 
will be referred to as the trapping region, see Fig. 1. 



Figure 1 
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4.2 Impact oscillator 

In the present Subsection we shall work in a general setting, proposed by Nord- 
mark |16j . The equation of motion is: 



d 2 x f dx 

dF = A [ x, dt' t 



(46) 



with periodic acceleration function: 



A{x >% t + 2n ) =A ( X >% t 



(47) 



The motion of the material point is confined to x > x c and the stationary 
impact boundary is given by x — x c . The impact law is of form = G (Uj), 
G (0) = 0, see for example Eq.( l3b[) . This equation of motion can describe linear 
impact oscillators as well as a ball bouncing on a vibrating table, see Section 
5.2.3 of [38] for the transformation relating systems with a moving boundary to 
systems with a stationary stop. 

The solution of (14U1) can be written in form: 



x = x (*,X c ,Uo,to) ' 



with initial conditions x (to, x c ,vo,to) = x c 



dx(t,x c ,vo,to) 



(48) 

= vq > 0, i.e. we 



assume that the impact occurred at t 
impact occurs at t = t\: 



to- Let us also assume that the next 



F (t!,x c , v ,t ) = x (ti,x c , v , t ) -x c = 0, 



(49) 



where the time of the second impact is defined by (|49j) as an implicit function 
t\ = f (x c , Do, to). Obviously, due to the initial condition, it follows that t\ = to 
is the trivial solution of (|49|) and can be thus ignored. We shall next assume 
that t\ — to is small - this assumption is clearly fulfilled in the phase of grazing. 
Expanding the function F we get: 



. dx 
^-^df, 



d 2 x 



+ O((i!-to) 3 )-^c-0. (50) 



Neglecting the trivial solution t\ = to we obtain from (|50l ) the approximate 
formula for the time interval elapsing between two subsequent impacts: 



to = 



-2v 
A ' 



(51) 



where Aq 



"IP, 



A (x c ,vo,tfi). 



10 



5 Basins of attraction 



We have studied dependence of the dynamics on the control parameters. The 
bifurcation diagram is shown in Fig. 2 for R = 0.85 and < A < 2.5. It can 
be seen that for A < X^} — 0.5094474 . . . only the grazing dynamics is possible. 

For A > X^} the first fixed point pj 1 ^ exists and is attractive and the second 

fixed point P s (2) exists and is attractive for A > A^ = 1.0188949 .... The first 
of 3 - cycles, which typically accompany 1 - cycles, appears for A = 1.78 while 
several small attractors can be seen for A > 1.5. 

Figure 2 

It follows that for A*/ < A < A^ there are two coexisting attractors only: 
the grazing manifold and the first fixed point P s . 

Figure 3 

Basins of attraction of these two attractors are shown in Fig. 3 above for R = 
0.85, A = 0.6. In the Figures 3 — 5 basins of grazing dynamics are yellow, basins 
of pj 1 ' are coloured red and basin of Pi 2 ' is green. Arrows denote direction 
of motion on stable and unstable manifolds of unstable fixed points P^p , P„ . 
Boundaries of the attractor's basins are determined by the nature of the unstable 
fixed point P^p . This unstable fixed point appears for growing A when A > 
X^r\ = ^jxj; ■ We have demonstrated in Section 3 that eigenvalues of the 
characteristic equation fulfill inequalities < Ai < 1 and 1 < A2 ^A > A^^J. 

Therefore the fixed point pfi is a saddle: it has one-dimensional unstable 
manifold W u 

(Pu^ as well as one-dimensional stable manifold W s 

(Pu ] ) in the 

phase space (A, V). The basin boundary is exactly the manifold W s (pfPj [55] . 

We can see in Fig. 3 fractalization of the basins for large velocities. This 
phenomenon is enhanced for larger values of A, see Fig. 4 below. 

Figure 4 

For higher values of A, A > A^ = — j-f]^-, the second pair of 1 - cycles, 

(2) (2) 

P s and Pu , exist (the first attractive and the second one a saddle). Basin 
boundaries are formed by W s an d W« ^P« 2 '^ , see Fig. 5. 

Figure 5 

In Figs. 3, 4, 5 fractalization of the basins of attraction is clearly seen. 
In other words, the manifolds W s (Pi 1} ),^(pi 2) ) , which contribute to the 
basins boundaries, become fragmented or shredded. The shredding process of 
the manifolds was predicted and qualitatively described by Whiston [T5] . It also 
seems that all three basins in Fig. 5 are intermingled in some parts of the phase 
space. 
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6 Discussion 



We have analyzed dynamics of the bouncing ball via the Poincare map. We 
have found approximate solution, (|39|) . of Eq. (|17aD defining time of the next 
impact as the implicit function, Xj+i = / (JQ, Vi), near the grazing manifold, 
except from some critical points, discussed below. It was possible to determine 
approximate form of the implicit function because the series expansion of (|17a[) 
converges near the grazing manifold, moreover the trivial solution, Xj+i = Xi, 
was recognized as an initial condition and was eliminated, accordingly. The 
solution (|39[) or (|40|) was shown be valid for more general impacting systems, see 
Eq. (|51[) . The critical points on the grazing manifold turn out to be zeroes of the 
post-impact acceleration. Let us note here that importance of the acceleration 
function was stressed by Whiston and Nordmark [TU [T51 [H] . 

We have also studied dependence on initial conditions for the map (TT7|) in 
the case of coexisting attractors, namely grazing dynamics and the first fixed 
point Ps , as well as the second fixed point P^ 2 \ It is interesting that for 
initial conditions with large velocities fractalization of basins boundaries can be 
observed. This effect was predicted by Whiston [15] . 

There are several possible directions for further investigations. Firstly, the 
role of critical points on the grazing manifold should be studied. Secondly, the 
shredding process, leading to fragmentation of basins of attraction, should be 
understood in detail. Finally, more general case of impacts can be considered. 
Namely, it is well known from experiments as well as from theoretical investiga- 
tions that the coefficient of restitution depends on the impact velocity, R = r (x) 
[40l 133] . There were several experimental as well as theoretical attempts to de- 
termine dependence of the coefficient of restitution on impact velocity, see for 
example [HI H2] . It is thus possible to use the approach described in the present 
work in the case of velocity dependent coefficient of restitution, substituting in 
the impact equation (|3b[) the constant R by a function r (i) . Alternatively, 
continuous model of low-speed collisions can be considered [33] . 
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Figure captions 



1. Fig. 1 

Trapping (medium solid line) and forbidden (dotted line) regions of the 
grazing manifold V = A cos X; V = j and V = sin X (thin solid lines), A = 
2.5. Critical points are at the ends of dotted lines. 

2. Fig. 2 

Bifurcation diagram, R = 0.85, < A < 2.5. Dotted lines enclose trian- 
gular region of grazing dynamics with —A < V < A. 

3. Fig. 3 

Basins of attraction, R = 0.85, A = 0.6. 

4. Fig. 4 

Basins of attraction, R = 0.85, A = 0.9. 

5. Fig. 5 

Basins of attraction, R = 0.85, A = 1.2. 
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